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We analyze the applicability of perturbative QCD (pQCD) approach to the is¬ 
sue of J/ip recombination at the Large Hadron Collider (LHC), and calculate the 
recombination cross section for cc recombination to form J/ip as a function of tem¬ 
perature. The charmonium wavefunction is obtained by employing a temperature 
dependent phenomenological potential between the cc pair. The temperature de¬ 
pendent formation time of charmonium is also employed in the current work. A set 
of coupled rate equations is established which incorporates color screening, gluonic 
dissociation, collisional damping and recombination of uncorrelated cc pair in the 
quark-gluon plasma (QGP) medium. The final J/ip suppression, thus determined as 
a function of centrality is compared with the ALICE experimental data at both mid 
and forward rapidity and CMS experimental data at mid rapidity obtained from the 
Large Hadron Collider (LHC) at center of mass energy yj snn = 2.76 TeV. 
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t"" Charmonium and bottomonium (quarkonium) 
]^iave been considered as important signatures to 
^|tudy the formation and properties of the quark 
Qrfuon plasma (QGP). Various experimental 
*Cind theoretical work j8-17| have been carried out 
• to study the charmonium and bottomonium sup- 
. press ion. These signatures are similar in many re¬ 
spects, but differ mainly in the aspect of secondary 
charmonium/bottomonium formation i.e. recombi¬ 
nation. Due to the heavy rest mass of the bottom 
quark, very little bottom quark and anti-quarks are 
formed even at the LHC energies. On the contrary, 
a relatively much larger number of charm quarks 
and anti-quarks are produced at the same energy, 
which can recombine to form charmonium even at 
the later stages of the QGP. There have been quite 
a few attempts for modeling the recombination phe¬ 
nomena to determine the effective quarkonium sup¬ 
pression. 


In Ref. 18 , the authors have used statistical 


hadronization to capture the J/ip recombination, 
while in 19, 20|, transport models have been uti- 
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lized. In the approach adopted in Ref. 21], re¬ 
combination cross section has been derived from the 
dissociation cross section using detailed balance ap¬ 
proach. In this work, we look at the applicability of 
pQCD to calculate the recombination cross section. 
While recombination is a non-perturbative process, 
we provide arguments and simulation results that 
show that pQCD can be a fair approximation at 
the LHC energies. We also model the temperature 
dependence of the recombination cross section. 


We model the recombination mechanism as a 
pair of uncorrelated charm quark and anti quark 
colliding to produce a gluon and a J/ip. The Feyn¬ 
man diagrams corresponding to this mechanism are 
given in Appendix A. To validate the applicability 
of pQCD, we investigate the momentum transfer in 
the collision in the cc center of mass frame of ref¬ 
erence. In the diagrams 1 to 5, it can be seen that 
the uncorrelated c and c pair annihilate to form a 
virtual gluon, which then creates a correlated cc 
pair and eventually forms a J/ip. The minimum 
energy of the incoming cc pair has to be equal to 
the J/ip rest mass = 3.1 GeV. In the center of mass 
frame, the net input momentum is 0, which implies 
that the virtuality of the gluon propagator needs to 
be > J/ip rest mass (3.1 GeV). This is within the 
realm of pQCD. In the diagrams 6 to 10, there is 
no annihilation of incoming cc pair, but there is a 
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momentum transfer via the virtual gluon. The mo¬ 
mentum of the bound quark and anti-quark in J/^ 
and the momentum of the outgoing gluon is trans¬ 
ferred via the gluon propagator. From the solution 
of the Schrodinger equation, it can be seen that the 
standard deviation of the momentum spread of the 
bound charm quark and anti-quark is about 0.5 GeV 
and this helps in reducing the error due to pQCD 
even if the emitted gluon is not very hard. The dis¬ 
tribution of the momentum of the outgoing gluon 
plays a significant role in determining the error due 
to pQCD. 

Given the large px range of many (though not all) 
of the outgoing J/i^ at the LHC, many of the gluons 
are likely be hard gluons. The accuracy of pQCD 
calculations would depend upon the percentage of 
hard gluons emitted. In the initial hard production, 
the correlated c and c are generated from the collid¬ 
ing nucleons, with the collision axis being along the 
beam axis. However, in the case of the formation 
of secondary J/i\) produced from recombination of 
the uncorrelated charm quark, and anti quark, the 
collision axis need not be directed along beam axis. 
Consequently, the px of the emerging J/-0 need not 
necessarily determine the hardness of the process. 
For instance, a pair of particles moving almost in 
parallel may undergo a soft collision, and could be 
Lorentz boosted due to their perpendicular (to the 
collision axis) momentum, resulting in a large px- 
But, since the collision axis and the J/ip emission 
can be directed in any random direction, the effect 
of Lorentz boost can be in either direction i.e., to 
either increase or decrease px■ Hence, though the 
px of the outgoing J/ip (gluon) in the lab. frame 
is not an exact measure of hardness, it can act as 
a probabilistic estimate for the momentum range of 
J/^ (gluon) in center of mass frame and thus can 
be an approximate measure of the hardness of the 
gluon emitted. In other words, a large momentum 
range in the lab. frame is likely to result in a large 
momentum in the center of mass frame. For Pb—Pb 
collision, at the ALICE mid rapidity, the px range 
is 0 < px < 8 GeV, while at the ALICE forward 
rapidity, the px range is px > 0. At the CMS mid 
rapidity, px range is 6.5 < px < 30 GeV. At forward 
rapidity, the longitudinal component of the momen¬ 
tum is significant, and can significantly increase the 
total magnitude of J/tf) momentum compared to 
ALICE mid rapidity. We describe this more quan¬ 
titatively in the latter part of the introduction. In 


the light of the foregoing discussion, the applicabil¬ 
ity of pQCD may increase in the order: 

ALICE (mid-rapidity) < ALICE(forward rapidity) 
< CMS (mid-rapidity,). 

In fact, from our simulation results, shown in section 
IV, it can be seen that the experimental J/ij) sup¬ 
pression data is reproduced to a much better extent 
for ALICE forward rapidity and CMS mid rapidity 
than at ALICE mid rapidity. 

The last vertex that needs to be explored is the 
outgoing gluon. The outgoing gluon can be a soft 
gluon, but as argued above, if the number of hard 
gluons are large, the discrepancy due to application 
of pQCD can be small. Secondly, in Appendix B, 
we show that the amplitude due to the soft gluon 
can partially cancel between the various pairs of di¬ 
agrams (Eg., between diagram 6 and 7 or 2 and 
3 etc.). If the soft gluon is perpendicular or par¬ 
allel to the collision axis, then it exactly cancels, 
and for other inclinations, it cancels partially to 
varying degrees. This, we expect, would soften the 
divergence of the Fermion propagator, and reduce 
the error due to divergence of the Fermion propa¬ 
gator in the pQCD calculations. To further vali¬ 
date all the above arguments for the approximate 
validity of pQCD, we did numerical simulation in 
which the amplitudes with gluon virtuality < 1 

GeV or outgoing gluon momentum < 1 GeV were 
zeroed out. The resulting recombination rate is in¬ 
dicated by the dashed line in Fig. 2. Comparing 
the difference between the solid and the correspond¬ 
ing dashed curves, the recombination decreases by 
a maximum of about 9% at ALICE mid rapidity 
and by a much smaller extent at ALICE forward 
rapidity. The difference is indicative of the number 
of soft gluons present. This corroborates with the 
above arguments on the gluons being mostly hard. 
To summarize, while an exact treatment of J/-0 re¬ 
combination would require a non-perturbative anal¬ 
ysis, the arguments and simulation results indicate 
sufficient grounds to believe that pQCD provides a 
fair approximation of the recombination mechanism 
at LHC energy. 

At ALICE forward rapidity, in order to estimate 
the longitudinal component, we analyze the total 
momentum for px > 0.35 GeV. At px = 0.35 GeV, 
the longitudinal momentum pi > 2.15 GeV for the 
rapidity range 2.5 < y < 4.0. The total J/^ mo¬ 
mentum becomes around 2.2 GeV, corresponding to 
energy = 3.8 GeV. The number of J/i^ particles less 
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than p T < 0.35 GeV may be small. As a support 


ing evidence, from Ref. 22], the J/-0 production 


peaks at around pt = 1.8 GeV at forward rapidity 
y/SistN = 7 TeV. Assuming the distribution to be 
similar for Pb — Pb collision at yTF/viv = 2.76 TeV, 
only a small number of JA0 particles will lie in the 
range pt < 0.35 GeV. Additionally, there would be 
J/-0 particles with pr < 0.35 GeV, but still with 
total momentum greater than 2.2 GeV. With the 
charm quark distribution (solid line) in Fig. 1, an 
estimate of the number of J/ij) particles with total 
momentum >2.2 GeV is 99%. Hence, for ALICE 
forward rapidity, we take 2.2 GeV as the lower mo¬ 
mentum cut off. 

In Ref. 


21 , the recombination cross section has 


been derived from the dissociation cross section 23 
using detailed balance approach. The dissocia¬ 
tion cross section developed in (23| is an extension 
over f24j by including confining string contribution 
using NRQCD. However, this formulation has been 
developed for soft gluons. One would thus expect 
it to be accurate for gluon momentum up to ~ 1.0 
GeV or 1.5 GeV. As a result, the resulting recombi¬ 
nation cross section would only be a rough approx¬ 
imation when applied to hard processes, especially 
at the ALICE forward rapidity (p > 2.2 GeV) and 
CMS mid rapidity {j>t > 6.5 GeV). The current 
pQCD formulation may be a better approximation 
at both ALICE forward rapidity and CMS mid ra¬ 
pidity. Even at ALICE mid-rapidity sufficient num- 



the application of temperature dependent formation 
time which is described later in section IIII Bl 


We use the pQCD calculations to determine the 
temperature dependent recombination time con¬ 
stant, r recom b{T) which is based on the temperature 
dependent cross section cr^ ecomb for cc —$■ J/-0 pro¬ 
cess at various energies of the cc pair. The charmo- 
nium wavefunction used in the calculation of crJ ecomb 
is determined by solving Schrodinger wave equation 
with a temperature dependent Debye color screened 
phenomenological potential. As mentioned earlier, 
we mainly consider the process of a quark and anti¬ 
quark annihilation to form J/ifj and a gluon in the 
final state, which is an a 3 process. In the a 2 process, 
where there is no final state gluon, the cc —>- J/-0 
reaction would involve only a very limited phase 


space domain, in which the sum of energy of in¬ 
coming cc pair needs to be exactly equal to the J/-0 
rest mass in the cc center of mass frame of reference. 
Due to this extremely limited phase space domain, 
it plays an insignificant role in the recombination 
process. The other a 3 processes involving gluon, 
quark and anti-quark in the initial state have simi¬ 
lar limited phase space restrictions. Moreover, they 
involve three body collision and these processes are 
ignored. The J/-0 dissociating into a cc pair, and 
cc recombining to produce J/ifj form a system of 
coupled rate equations, with each feeding into the 
other. In this work, we model this phenomenon by 
a system of coupled rate equations. 

Furthermore, Jand cc would be formed from 
the interaction of light quarks and gluons in the 
medium. The pQCD framework developed here has 
used a relativistic modeling of the initial cc pair, 
and hence can be utilized to calculate the temper¬ 
ature dependent cross section for the lighter quarks 
to form J/'ijj. In section II (D), we analyze the im¬ 
pact of the light quarks. Additionally, some of the 
cc would also decay back to the lighter quarks and 
gluons. We ignore this decay process in the current 
work. It is to be noted that the process involving 
the lighter quarks, would not be part of the feed¬ 
back mechanism happening between J/-0 and cc, 
and hence, they would have a much lesser impact 
to the overall suppression or enhancement (via re¬ 
combination) as compared to the cc -H- J/ip process. 
We also ignore the running of the coupling constant, 
and instead use a fixed value of the coupling con¬ 
stant. 

It is found that the charmonium suppression in 
the QGP is not the consequence of a single mecha¬ 
nism, but is a complex interplay of multiple mecha¬ 
nisms like color screening, gluonic dissociation, col- 
lisional damping and recombination. We develop 
a framework of rate equations, which incorporates 
all the above mechanisms. Color screening is based 
on quasi-particle model (QPM) equation of state 
(EOS) of the QGP described in j25j. The QGP 
is expanding under Bjorken’s scaling law, applica¬ 
ble mainly at mid rapidity. Due to the modifica¬ 
tion of the formation time by temperature j26[, we 
show that the color screening becomes very small. 
Hence, we ignore color screening to compute sup¬ 
pression at forward rapidity. This could lead to 
a little under suppression in peripheral collisions. 
The process of gluonic dissociation and collisional 
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damping is based on the formulation developed by 
Wolschin j23[. At the LHC energies, absorption is 
expected to be negligible, since the cc pair would 
behave almost as a color singlet while traversing 
through the nucleus and hence have negligible in¬ 
teraction with the nucleus. The p t broadening due 
to Cronin effect is inconsequential while explaining 
the p t integrated suppression at the mid rapidity 
for various centrality bins. That is why, we incor¬ 
porate only shadowing based on the work done by 
Vogt [27[ as a CNM effect in our present calcula¬ 
tion. In fact, in Ref. [27[, the shadowing for Pb—Pb 
collisions at LHC energy y/s^N = 5.5 TeV was de¬ 
termined. We use the same framework to calculate 
shadowing at y /snn = 2.76 TeV energy. Our results 
for J/-0 suppression versus N part are compared with 
the experimental data in the forward and mid rapid¬ 
ity region obtained from the ALICE (d, 5j and mid 
rapidity data from CMS experiments at LHC 
energy y/s^N = 2.76 TeV. We find that our predic¬ 
tions show a reasonably good agreement with the 
suppression data, particularly at ALICE forward ra¬ 
pidity and CMS mid rapidity. 

The organization of the rest of the paper is as 
follows. Section II introduces the system of coupled 
rate equations, and describes the pQCD calculation 
of the temperature dependent recombination cross 
section. Section III briefly describes the mechanism 
of color screening, gluonic dissociation, collisional 
damping and CNM effects. The last part of section 
III, describes the inclusion of all these mechanisms 
in the rate equations. The section IV gives the re¬ 
sults and discussions and finally conclusion is given 
in section V. 


II. RATE EQUATIONS AND THE 
RECOMBINATION CROSS SECTION 

A. The coupled rate equations 

Once the charmonium bound state is formed af¬ 
ter the formation time, we assume that there are 
two reversible processes which are in play. The first 
one is the dissociation of charmonium into its con¬ 
stituent charm and anti-charm quark. The second 
is the recombination of the charm and anti-charm 
quark to again form charmonium. For a QGP with 
instantaneous volume V (t), we model these two pro¬ 


cesses using the following set of rate equations: 


dN JH) = _ Tdiss{T)Nj/ ^ t) + r reeomfe (r) iVc( T t / ) ^ 5 (t) ( i) 


dt 


dN c 

dt 


— r diss(T)Nj/jp(t) — r recomb(T ) 


V(t) 
N c (t)N- c (t ) 
V{t) 


where R^s and T recom b are J/%b dissociation and cc 
recombination rate, respectively. The above equa¬ 
tions are solved numerically. The initial conditions 
are given by Nj/^( 0), A r c (0) and N s ( 0) after includ¬ 
ing shadowing. N c ( 0) and W(0) are determined 


from 


and taken as 75 based on the MRST HO 
structure function of nucleon for the most central 
collision. The ratio of Nj P „,. to N pp is determined 
from the experimental data 


29 . The value of N pp 


and N pp is obtained from the open charm D meson 
data from LHC j29| and extrapolated till p t = 30 


GeV using the techniques mentioned in |30|, 131 


The open charm D meson distribution is then scaled 
along the momentum axis according to the momen¬ 
tum fraction of the charm quark present in the D 
meson to obtain the charm quark momentum dis¬ 
tribution. The resultant normalized charm distri¬ 
bution is shown in Fig. 1 (solid line). The extrapo¬ 
lated value of N C1 after adjusting for the difference 
in luminosity for N pp ^ and N pp data, is then used 


n pp 

to determine the ratio -j$r- This ratio is then used 
to obtain AGm( 0) from iV c (0). The obtained val¬ 
ues of W(0) and Nj/^(0) are then scaled for various 
centrality bins according to Taa and shadowing to 
arrive at the initial conditions for the rate equa¬ 
tions. Modeling of the shadowing is described in 
the section on CNM effects. 

In view of the fact that the initial number of J/-0 
is negligible, the coupled rate equations can be sim¬ 
plified leading to an analytical solution as outlined 
in [2lj]. The difference between the numerical solu¬ 
tion of the coupled rate equation and the analytical 
solution of the simplified equations is small. In this 
work, we have used the coupled rate equations for 
its more generic and allows the analysis of other ef¬ 
fects, for example, the effect of light quarks in sec¬ 
tion QT]21 It also helps in more accurately analyzing 
the discrepancy in J/ip yield due to uncertainty in 
the initial conditions as discussed in Section m 

The value of recombination time constant (re¬ 
combination reactivity) is given by T recomb {T) = 


f dp 1 dp 2 f c (pi) fc(p2)Vrecomb(Pl,P2)v re l, where V rd 
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is the relative velocity between c and c in the lab¬ 
oratory frame. f c (Pi) and /c(p 2 ) are distribution 
function of c and c with four momenta p\ and p 2 , 
respectively. crJ' ecomb is computed using pQCD cal¬ 
culations. The suppression or enhancement mod¬ 
eled by the rate equations is defined by 


Nj/^(t QG p) 


( 2 ) 


where tQGP is the lifetime of QGP (= 5 fin) at the 
LffC energy, and to is some initial time, which is 
taken to be the thermalization time = 0.5 fm. 


wavefunction ipT^k, k'). The charmonium bound 
state is modeled non-relativistically (a solution of 
a Schrodinger equation) [2ll . 33}, with the charm 
quark and anti-charm quark assumed to have very 
small momentum spread. From our solution to the 
Schrodinger equation at the minimum QGP tem¬ 
perature of 170 MeV, i.e. before hadronization, the 
J/ip wavefunction is seen to have a one standard 
deviation spread of 0.5 GeV (approx.). At higher 
temperatures, the wavefunction expands spatially 
leading to a decrease in the momentum spread. 

The four momentum conserving delta function 
5 4 (pi + P 2 — k g — k — k’) would then imply that 
the momentum spread of bound state charm quarks 
and anti-charm quarks needs to be compensated by 
the momentum spread of input quarks and outgoing 
gluon. The incoming quarks and anti-quarks may 
naturally acquire a momenta spread due to confine¬ 
ment in QGP medium. The remaining needs to 
be compensated by the gluon, which may make the 
gluon go off shell by a small amount (<0.5 GeV). 
As an approximation, we ignore the effect of the 
momentum spread of the incoming charm quarks, 
anti-charm quarks and the outgoing gluon and place 
the gluon fully on shell. This could however be a 
source of error for soft gluons emitted particularly 
at ALICE mid rapidity. 

With the above approximation of the input quark 
being in a definite momentum eigenstate and the 
gluon being on shell, the probability amplitude for 
the various Feynman diagrams (given in Appendix 
A) is given by 


B. pQCD calculation of (J recom b 

We use pQCD to compute the cross section of 
charm quark and anti-charm quark to form a fi¬ 
nal state gluon and J/i^. In Ref. [32}, Feynman 
diagrams for QQ scattering to order a 2 , and for 
QQ annihilation (production) to (from) light quarks 
and gluons are given. Since recombination involves 
production of J/-0 from uncorrelated cc pair (in¬ 
stead of light quarks and gluons), the Feynman di¬ 
agrams in [32[ need to be modified. The Feynman 
diagrams corresponding to the recombination pro¬ 
cess are given in the Appendix A. We treat the J/^j 
bound state as a linear superposition of definite mo¬ 
mentum eigenstates of a free particle. The contri¬ 
bution of each momentum eigenstate is given by the 


1. Mi = (v(p 2 )ig'Y„t a u(pi)^( 1 ^.gf abc ) 

g v »(k l +k g y+g^(-k g -k 2 ) v +g™(k2-k 1 )i*) 

j%tfi(k g ) (/ u(k){ig^ a t c )v{k')^ T (k,k')dk^j 

2. M 2 = (v(p 2 )(ig^X)y^(ig^t b )e^\k g )u(p 1 )^ 
(nip) (/ u(k)(ig-i a t a )v{k')^ T {k,k')dk^ 

3. M 3 = ( y(p 2 )(ig^^t b )e^ b (k g )^ i (ig^t a )u(p 1 )) 
(nip) (j u(k)(igj a t a )v(k')'il; T (k,k')dk\ 

4. M 4 = (u(p 2 )(^7^“MPi))(^fp) 
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( J u(k)(ig^^t b ) m (ig'yX)v(tf) 

tf> T (k,k')dlc)e ,l ’ b (kg) 

5. M 5 = (v(p 2 )(ig^t a )u(p 1 j S j(^^-) 

( J u(k){ig^ a t a ) ^ l _ m {ig~i^t b )v{k') 

‘th(k,k , )dk > jef t ’ b (k g ) 

6 . M 6 — J f ^u(p 2 ) {igivt a )v{k '){—^ 

(u(k)(ig-i a t a ) ^ m (^7//MPi)^) 
rj} T {K k f ) j dke^ikg) 

7. M 7 = j (v{p 2 ){igivf) ^ vt a )v{k')^ 

( ^ ) (u(k)(ig^X)u(piX T (k, fc')) dk e^ b (k g 


is not possible. The above amplitudes are evalu¬ 
ated explicitly. The averaged M 2 is then obtained 
as Haii spins [(S-£i M*) (E \t\Mi)]. The factor ^ 
comes from averaging over the spins and colors of 
the incoming charm quark and anti-charm quark. 
In order to avoid the ghost terms, only the trans¬ 
verse component of the gluon polarization term e M,fc 
is used. 

The charmonium wavefunction ipx(k , k r ) at tem¬ 
perature T is obtained by solving the Schrodinger 
equation with a phenomenological potential using 
the real part of the potential given in Eq. Q] m c 
is the mass of charm quark, while u and v are the 
respective spinors of c and c quarks, respectively. 
The value of charm quark mass, m c , used in the 
pQCD calculations as well as in other calculations 
(e.g., potential model computations) in this paper 

is 1.275 GeV. We note that k = ( \J\k\ 2 + m 2 , k) 

and k! = (\f \k\ 2 + m 2 , — k ) as 4-momentum of the 
charm and anti-charm quarks. The recombination 
cross section in the laboratory frame is then com¬ 
puted in the cylindrical coordinates. 


- 2 tt 


r*2-rr 




8 . M 8 — / v(p 2 ){ig'iX)v{k')( 


P{PhP2)-V re l = ~ 


do 


d/3- 


kl 


^ ^ rlh com 

~ aK a\\ 


u(k)(ig^t b ) — 

y rfl — Tn 

Mk,k') dke^ b (k g ) 


-(igiX)u(pi) 


2 vr J o .A, 4vr 2 . 

M 2 (pi,p 2 , kg, k, kj/^j Q(k l j /r k h j^) k g± 


x 


2 E 1 2 E 2 Ej^ E g h lab 



(3) 




9. Mo = 


v(P2){igiX)- 


— m 


{Wlpt )v(k') 


u(k)(ig^X)u(piXT{k,k')) )dk 


**'(kg.) 


10 . Mio — 



v(p 2 )(igiX)v{k')gf bac 


g^i-kg-kX+g^X-kX+g^X+kg ) 1 
1 


For ALICE suppression data at forward rapidity, 
we take k l j^ to be 2.2 GeV and fcj/q to be a very 
large number, while at ALICE mid rapidity, the 
values are 0 and 8 GeV. For the CMS suppres¬ 
sion data, the values of fcW and k b ^ are 6.5 and 
30 GeV, respectively. The cylindrical coordinates 
are chosen so that k g jj m is along the direction of 
velocity of the center-of-mass ( jj com ) of the cc sys¬ 
tem. The angle (3 in the cylindrical coordinates is 
given by the angle of rotation around the v com axis, 
and 7 = 


h 2 h 2 ' 

r» 2 -r»l 


uXiig^XuXXrik, k')) dk e^ b (k g 


Here g is the strong coupling constant. The 
wavefunction X{k, k') is obtained numerically, and 
hence analytical evaluation of the above expressions 


—r — -. k a , kj m, pi and p 2 are con- 

j strained by the 4-momentum conserving delta func¬ 
tion <5 4 (pi + p 2 — k g — kj/^). In general, for all the 
quantities, we use the notation k and p to denote 
the spatial component of the 4-momentum k and p. 
h lab is the outcome of integrating the energy con¬ 
serving delta function S(Ei + E 2 — E g — Ej/^), and 























7 


TABLE I: Temperature dependent values of Y recomb 


Temperature 

(MeV) 

Pt <8 GeV 
(mb) 

p > 2.2 GeV 
(mb) 

6.5 < pt 
(mb) 

0 

0.0000 

0.0000 

0.0000 

50 

0.3624 

0.0240 

0.0000 

100 

0.5697 

0.1275 

0.0003 

150 

0.4823 

0.1698 

0.0017 

200 

0.3706 

0.1649 

0.0040 

250 

0.2378 

0.1222 

0.0052 

300 

0.1219 

0.0692 

0.0044 

350 

0.0389 

0.0238 

0.0020 

368 

0.0224 

0.0140 

0.0013 


is given by 

] lab _ i J_ Eg _ v corn cos( C) 

L 1 ' (■yEy’™.(l+v com )' 2 ) i + ^com 

( is the angle between kj™ and v com and 6 is the 
angle between c and c 3-momenta in the laboratory 
frame. The superscript ” com” refers to the input cc 
center-of-mass frame. 

2 

We take the value of a = = 0.31. The value 

of T recomb obtained as a function of temperature is 
shown in Fig. 2. The numerical values correspond¬ 
ing to Fig. 2 is given in Table I. The value of T recom b 
at around 200 MeV (180 MeV to 200 MeV is about 
the lowest temperature of the QGP, and the recom¬ 
bination is likely to be highest at this temperature) 
is about 0.37 mb at ALICE mid rapidity, 0.16 mb 
at ALICE forward rapidity and 0.004 mb at CMS. 
The range 0.48 mb to 0.37 mb between the tempera¬ 
tures 150 MeV to 200 MeV, at ALICE mid rapidity 
is comparable with T recomb = 0.65 mb used in [301 ]. 
but on the lower side. The dashed lines indicate 
the value of T recomb , when the amplitudes are ze¬ 
roed out for gluon propagator virtuality or outgoing 
gluon momentum < 1 GeV. The difference between 
the solid line and the corresponding dashed line is 
indicative of the number of soft gluons present. 

Figure 3 shows the average value of a recomb .v re i 
for the special case of |pi| = |p 2 | at 200 MeV. Apart 
from the initial increase, the value of a recomb de¬ 
creases with the center-of-mass energy. The a recomb 
for the ALICE mid rapidity is higher than the value 
for ALICE forward rapidity which is higher than the 
value for CMS mid rapidity. 

The decreasing value with energy is consistent 
with the fact that recombination decreases with in¬ 
creasing energy, as can be seen from the experimcn- 



FIG. 2: Y recomb as a function of temperature. 

The legend is: 

A = T recomb @ ALICE mid rapidity; 

B = T recomb @ ALICE forward rapidity; 

C r recomb @ CMS; 

D r recomb @ ALICE mid rapidity, e = 1 GeV] 

E = Y recomb @ ALICE forward rapidity, e = 1 GeV 
e = momentum cut off for virtual gluon 
propagator and outgoing gluon. 

tal data in Figs. 11, 12 and 13, where the observed 
recombination is highest at ALICE mid rapidity 
(lowest momentum range). 

1. Factorizability 

The formation of J/ip from cc pair is a non- 
perturbative effect involving exchange of soft glu¬ 
ons and occurs over a long time scale. Referring to 
Eq. B17 in Appendix B, for soft outgoing gluons, 
the infrared cancellation is only partial. If p\ or p 2 
(the spatial components of pi and p 2 respectively) 
were zero, then the cancellation would have been 
exact. For the case of soft outgoing gluons, p\ or 
P 2 is small but not zero. Hence, the cancellation of 
infrared divergence is partial. Consequently for soft 
gluons, factorization between the potential model 
based J/if; wavefunction and cc scattering would 
lead to inaccuracy in the result. For very hard out¬ 
going gluons, for the cc scattering sub-process, one 
can ignore the momentum spread of the J/-0 wave- 
function which is much smaller in comparison to 




















Total center of mass Energy (GeV) 


FIG. 3: r rec0 mb as a function of total energy of cc 
at 200 MeV 


the gluon momentum (J/ip momentum spread ~ 
0.5 GeV at temperature of 170 MeV). This means 
that the scattering of the c and c occurs over a rela¬ 
tively short time scale. Hence, for very hard gluons, 
it is possible to factorize the potential model based 
J/rj) wavefunction and the hard cc scattering pro¬ 
cess. At least, to the order a 3 in the perturbation 
theory, it would then be possible to represent: 
cr(cc —* J/if> + g) —* cr(cc —* cc + <?)^ 2 (0), 
where "0(0) is the value of the J/ifj wavefunction at 
0 spatial separation and captures the long duration 
interactions. In the light of the above arguments, 
the validity of factorization of the recombination 
cross section would decrease in the order 
CMS mid rapidity (p T > 6.5 GeV) > ALICE for¬ 
ward rapidity (p > 2.2 GeV) > ALICE mid rapidity 
(p T > 0 GeV). 

In this work, in order to maximize the accuracy at 
ALICE conditions, where the recombination is high, 
factorization is not employed to determine the re¬ 
combination cross section. 


C. Volume scaling computation 


The volume scaling is based on the quasi-particle 
model (QPM) equation of state (EOS) of the QGP 
and the concept of constant entropy condition. The 
QPM EOS used for computing the volume scaling 


is described in 25 


In the QPM, the Reynolds number is given by 
R^ 1 = (4.0 rj)/ (3.0 T 0 7 o s); where rj = j-, entropy 
density s = 16.41 GeV 3 and T 0 is the temperature 
at t 0 = 0.5 fm. The volume is then given as: 


V(t) = v 0 (t 0 /t ) 


R- 1 - 1. 


where Vq is the initial volume at r 0 , and is taken as 
v 0 = 7 rRp b r 0 , with R Pl , = 6.62 fm from the Woods- 
Saxon distribution 34 . 


D. Lighter quarks and gluons 


The lighter quarks and gluons would also form 
J/-0 and cc. After including these processes, the 
rate equations get modified as 


dNj/^ 

dt 

dNc 

dt 


^'diss^J/ipit') T r recomb 


^ diss-^j/i/j (t) r 

recomb 


N c (t)N s (t) 

V(t) 

N c (t)N- c (t) 

V(t) 

Rq,g^rcci 


where 


A q,g^J/p = I dp\dp2 f p (P\)fp(;P2) X 

lg^J/^{Pl,P2)vZ(Pl,P2) 


cr 


A 


q,g->cc 


5 = / dp x dp 2 


fp{Pl)fp(P2) X 


cr 


q,9->cc 


(Pl,P2)v%l(Pl,P2) 


The light quarks and gluons are assumed to 
be thermalized and hence momentum distribution 
fpi.Pi) an d fp{p2) are taken to be Fermi Dirac 
distribution . . ; )+1 for light quarks and anti¬ 
quarks and Bose Einstein - . i/T) _ [ for the glu¬ 

ons. The baryonic chemical potential p is taken to 
be vanishingly small at the LHC. The value of de¬ 
generacy factor for light quarks (g q ) and for gluons 
(g g ) is taken as 6 and 16, respectively. 

For the lighter quarks, only the first five of the 
ten Feynman diagrams are applicable. With the 
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first five diagrams, the value of & g P 2 ) for 

the lighter quarks are not expected to be negligible 
when compared to that of cc —* J /-0 cross section. 

However, with vanishing value of /i, it is seen 
that most of the quarks and gluons lie in the low 
energy region. The semi-log plot in Fig. 4 shows 
that at around 1.5 GeV, the particle density reduces 
to around 10 -5 GeV -4 . This is not sufficient to form 



FIG. 4: Fermi Dirac and Bose Einstein distribution 


either J /-0 or cc in any significant quantity, resulting 
in a negligible value for A q^j/^p and A qtg ^ cE . In 
view of the above facts, with a Fermi-Dirac or Bose- 
Einstein distribution for a thermalized quarks and 
gluon, we neglect the contribution arising due to the 
light quarks and gluons to the rate equation. 


III. GLUONIC-DISSOCIATION, 
COLLISIONAL DAMPING, CNM AND 
COLOR SCREENING 


A. Gluonic-dissociation, Collisional damping 
and CNM effect 

The modeling of Gluonic-dissociation, collisional 
damping are on the same lines as [2l[. For ease of 
reference, we reproduce the quark, antiquark singlet 


potential used in this work 23 . 


V (r,m D ) = ?^(l-e- m ° r ) 


m D 


D -m n r 


a eff 1 1Jl D + 


( 4 ) 


i a effT 


r°° 2 zdz 

L (i +* 2 ) 2 


sin (?«£> r z ) 
rri£,r z 


where mo is the Debye mass given by = 

Nf = 3 = number of flavors; 


Ana? ( yy + qy 


aj = 0.49; a string = 0.192 GeV 2 . a e ff = -y, where 
we have taken a = 0.22. This value of alpha, along 
with aj = 0.49, gives the dissociation temperatures 
close to the values 381 MeV for J/-0, 190 MeV for 


i/j', and 197 MeV for Xc 35 


The CNM effect is also on the same lines as 21 


The CNM formulation in the current work and in 
Ref. 


is based on the work done by 


36 . The 


rapidity distribution for the computation of X\ and 
X 2 is taken to be \y\ < 0.9 for ALICE mid rapid¬ 
ity, 2.5 < \y\ < A for ALICE forward rapidity and 
\y\ < 2.4 for CMS. The shadowing parametrization 
EPS09 [37] is used to obtain the shadowing func¬ 
tion, and the The gluon distribution function in a 
proton has been estimated using CTEQ 6 |3£ 


B. Color Screening 

The_color screening model in this work differs 
from [ 2 l[ in terms of the use of temperature de¬ 
pendent formation time. The color screening model 
used in the present work and in Ref. j2l[ is based 
on pressure profile [35|] in the transverse plane and 
cooling law for pressure based on QPM EOS [H, 33! 
for QGP. The cooling law for pressure is given by 


, . A B C D 

p r,r = A+ — + - + —a, 

Tl t t c » 


( 5 ) 


where A = — ci, B = c 2 c 2 , C = 3 ^ 2 !^ and D = c 3 . 
The constants ci, c 2 , c 3 are given by 

_ „ _ ,, 4?? . 

• ci - -c 2 r - 3 ^ 7 , 


c 2 = 


£ 477/I 1 


c 3 

4r? 

T 


cl- 1 


T a 


C Po + Cl) Tq‘ 
(^- 1 ) 


2 —1 
c 2 c 2 r 0 
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The above constants are determined by using dif¬ 
ferent boundary conditions on pressure and energy 


density described in 251, |33 


Writing the above equations at initial time r = 
To and screening time r = r s and combining with 
the pressure profile [25], we get the following two 
equations: 


B C 


D 


p(r 0 , r) = A + — + — + —j= p(t 0 , 0) h(r); ( 6 ) 


r. 


T 0 


BCD 

p(Ts , r) = A + — + — + —2 = Pqgp, (7) 

Ts T s T 

where Pqgp is the pressure of QGP inside the 
screening region required to dissociate a particu¬ 
lar charmonium state and it is determined by QPM 
EOS for QGP medium. The above equations are 
solved numerically and the screening time t s is 
equated with the charmonium formation time tf 
= 'yTf(T) to determine radius of screening region 


r, [25|, |26|, [33]. 7 is the Lorentz factor corresponding 


to the transverse energy of J/ip. The temperature 
dependent formation time fp(T) is determined us¬ 
ing the simulation of time independent Schrodinger 
equation as [ 2 f 


Tf{T 1 ) = Tf(T 0 ) 


Ebind(T 0 ) 
Ebind(Tl ) 


( 8 ) 


Figure 5 depicts the formation time for J/-0 ob¬ 
tained from both the time independent Schrodinger 
equation and the time dependent Schrodinger equa¬ 
tion. The temperature associated with the initial 


conditions, T init 26], is taken as a value between 200 


MeV and 150 MeV based on the centrality. The ex¬ 
pression for survival probability due to color screen¬ 
ing can be obtained as: 


S c (pt , N part ) 


2(a + 1 ) [ Rt 


tt R 2 t 


dr r (pmax (r) 


B 2 r 


where a = 0.5, Rt and (which is a function of 


p t and r s ) are defined in 25|, [33 


C. Final suppression 


In the computation of the final suppression, all 
the above mentioned effects are included. 



FIG. 5: Comparison of the J/X> formation time 
between the time dependent Schrodinger equation 
(TDS), the time independent Schrodinger equation 
(TIS) and the vacuum formation time. 


In Fig. 2, it is shown that the recombination rate 
r recomb decreases with temperature at high temper¬ 
ature range, again as a consequence of Debye mass 
rriD widening (spatially) the J/ip waveform. This 
raises the question whether the decreased recom¬ 
bination rate is remodeling the suppression due to 
color screening mechanism described in section III 
(D). The temperature dependent color screened re¬ 
combination cross section is based on the interaction 
of uncorrelated charm quark and anti-charm quark 
produced in the latter stages of the QGP, and thus is 
unrelated to the color screening mechanism (as de¬ 
scribed in section III (D) ), which acts upon the cor¬ 
related charm and anti-charm quark pair produced 
during the initial hard collisions. This implies that 
the processes of color screening and recombination 
apply to the different non-overlapping population of 
charm quarks (correlated and uncorrelated). Hence, 
both mechanisms are required to model the charmo¬ 
nium suppression in the QGP. 

The effect of shadowing and the rate equation 
involving gluonic dissociation, collisional damp¬ 
ing and recombination, can be put together as: 
ScnmSdr • The question arises as to how the suppres¬ 
sion due to the color screening can be combined. 
Given the large temperature dependent formation 
time, the process of gluonic dissociation and colli- 
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sional damping can happen even before J/rj} is fully 
formed. Hence, the processes of color screening and 
gluonic dissociation and collisional damping may 
overlap in time. 

In order to incorporate the color screening into 
the system of rate equations within the framework 
of the above arguments, we first determine the 
equivalent T c corresponding to S c . 


r c — — \n(S c ) / (tQGP — t 0 ). 


We then have T diss = T dg +T c , where T dg = T 


dg 


diss.nl 


+ 


r damp, with r disS) ni and T damp being gluodissociation 
and collisional damping dissociation rates jUl, 26 


At ALICE forward rapidity conditions, due to the 
non-applicability of Bjorken’s scaling law at forward 
rapidity, we ignore the color screening mechanism. 


at the ALICE mid rapidity experiment, which has 
the lowest momentum range. 

The predicted J/ij) suppression due to color 
screening alone is shown in Figs. 6 and 7. The 
experimental data on J/'ip suppression for ALICE 
mid rapidity [5j and CMS j(3] at the LHC energy 
\Jsnn = 2.76 TeV versus N part is also shown for 
comparison. Similar to the case of T, with temper¬ 


ature dependent formation time 26], the effect of 


suppression due to color screening mechanism be¬ 
comes very small or almost disappears. This result 
forms the basis of ignoring color screening at AL¬ 
ICE forward rapidity. Color screening with either 
temperature dependent formation time or constant 
formation time does not capture the suppression for 
both ALICE and CMS data simultaneously. This 
indicates some other processes may be playing a 
role. 


IV. RESULTS AND DISCUSSIONS 


The experimental charm quark distribution in 
Fig. 1 is a static distribution, and does not 
vary with temperature or time. To overcome this 
difficulty, the modified Fermi Dirac distribution 
N (X) / (exp(XE / kT ) + 1) is used, where N(X) is a 
normalization factor. The value of A is fixed by com¬ 
paring the distribution to the experimentally deter¬ 
mined distribution at temperature T c = 180 MeV. 
The modified Fermi Dirac distribution is employed 
in all the simulations. The Tsallis distribution, 


N(X, q) [1 + (1 — q)(—XE/T)}~ , with q -> 1+, 
gives very similar results (not shown in this paper). 
Our calculation of the temperature dependent re¬ 
combination cross-section is shown in Fig. 2. As the 
temperature approaches the dissociation tempera¬ 
ture of J/-0, the value of T recom b approaches zero in 
accordance with the expectation. At very low tem¬ 
perature, near 0 MeV, the charm quark and anti¬ 
quark do not have sufficient energy to form J/-0, 
and hence, the V recomb approaches 0 at low tempera¬ 
ture also. Qualitatively, the temperature behavior is 
similar to the temperature dependence obtained in 
Ref. [39] for be, be bound states. The relevant tem¬ 
perature for QGP would be above T c = 180 MeV. 
It is also seen from Fig. 3, that at high momentum 
range at the CMS, the cross section is very small 
compared to the cross section at the lower momen¬ 
tum range of the ALICE experiments at both for¬ 
ward and mid rapidity. The cross section is highest 



FIG. 6: Our predicted suppression due to only 
color screening. Experimental data on J/-0 
suppression at ALICE (mid rapidity) LHC are 
taken from j&j. 


Figure 8 shows the effect of shadowing, dissocia¬ 
tion and recombination for the ALICE mid rapidity 
(low momentum range). The curve ’’only dissoci¬ 
ation” is obtained by assigning Y recornb = 0. The 
difference between the ’’only dissociation” and ’’re¬ 
combination and dissociation” curves is the result 
of recombination. 

It can be seen that the recombination increases 
with centrality. Recombination depends on the 
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FIG. 7: Our predicted suppression due to only 
color screening. Experimental data on JM 
suppression at CMS LHC are taken from [6 ; ]. 


FIG. 9: Comparison between the effect of color 
screening, recombination and shadowing at 
momentum values corresponding to ALICE data 
at forward rapidity. 



FIG. 8: Comparison between the effect of color 
screening, recombination and shadowing at 
momentum values corresponding to ALICE data 
at mid rapidity. 


product N c N 5 . If we set N c = N 5 , recombination 
depends quadratically on N c . On the other hand 
dissociation depends linearly on Nj/^. In the pe¬ 
ripheral collisions both Nj and N c are small, but 
due to the quadratic dependence on N c recombina¬ 
tion is lower in the peripheral collisions. Figures 9 



FIG. 10: Comparison between the effect of color 
screening, recombination and shadowing at 
momentum values corresponding to CMS data. 
The curves ’’recombination and dissociation” and 
’’only dissociation” almost overlap, indicating 
negligible recombination. 


and 10 show the effect of shadowing, dissociation 
and recombination at ALICE at forward rapidity 
and at CMS at mid rapidity. Recombination is seen 
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to be smaller at ALICE forward rapidity than at 
ALICE mid rapidity, and becomes negligible at the 
CMS where the momentum range is very high. The 
dependence of recombination with centrality at AL¬ 
ICE forward rapidity remains similar to ALICE mid 
rapidity. 

We now compare the final suppression data of 
J/^ at ALICE mid and forward rapidity in Figs. 
11 and 12. In the central region where recombina¬ 
tion is highest, it can be seen that J/ip suppression 
at forward rapidity is in reasonable agreement with 
experimental data, but at ALICE mid rapidity, the 
suppression is overestimated due to the recombina¬ 
tion being underestimated. As discussed earlier, the 
longitudinal component of momentum gives rise to 
higher momentum at forward rapidity as compared 
to mid rapidity for the same pr- Due to the larger 
momentum, pQCD calculations at forward rapid¬ 
ity are more reliable than at mid rapidity. Clearly 
non-perturbative QCD is required to explain the re¬ 
combination data more accurately at ALICE mid¬ 
rapidity. 

Finally, Fig. 13 shows the final suppression of 
J/"0, which is comparable to the experimental CMS 
data [bj up to a large extent throughout the whole 
range of N part . Recombination was seen to be al¬ 
most vanishing for the CMS momentum range in 
Fig. 10. 

In summary, we have been able to explain to a 
reasonable good extent, the J/-0 suppression at AL¬ 
ICE forward rapidity and CMS mid rapidity. Even 
though the Jj\b suppression data has not been ex¬ 
plained at ALICE mid rapidity accurately, it can 
be inferred from Figs. 8, 9 and 10, that pQCD 
correctly predicts the decreasing trend in recombi¬ 
nation from ALICE mid rapidity to ALICE forward 
rapidity to CMS mid rapidity. 

Comparison of Fig. [13] with Fig. 10 in Ref. [2l| 
shows similar suppression of J/-0. However, the fig¬ 
ure in Ref. j2lj indicates that the recombination is 
more than 20% for certain centrality bins. The re¬ 
combination cross section based on pQCD calcula¬ 
tion results in only about 2% increase in relative 
J/^ yield due to recombination at CMS mid ra¬ 
pidity (Fig. [TO]). Experimental data relating J/i[) 
Raa to pr at the CMS [40] indicate recombination 
may become very small at high px- The increase 
in suppression due to lesser recombination is offset 
by a decrease in suppression due to color screen¬ 
ing due based on temperature dependent formation 



FIG. 11: Final suppression with color screening, 
gluonic dissociation along with collisional 
damping, recombination and shadowing. 
Experimental data on suppression at ALICE at 
mid rapidity are taken from (5j. 



FIG. 12: Final suppression with color screening, 
gluonic dissociation along with collisional 
damping, recombination and shadowing. 
Experimental data on suppression at ALICE at 
forward rapidity are taken from [3]. 

time. Coincidentally, the net J/^ suppression at 
CMS mid rapidity is similar between the current 
work and |21h 
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FIG. 13: Final suppression with color screening, 
gluonic dissociation along with collisional damping, 
recombination and shadowing. Experimental data 
on suppression at CMS are taken from [Gj]. 


Apart from the discrepancy due to pQCD ap¬ 
proximation, an additional source of uncertainty 
in suppression might be that strange quarks may 
not be fully thermalized. If we compare the non- 
thermalized distribution in Fig. 1 and the thermal¬ 
ized Fermi Dirac or Bose Einstein distribution in 
Fig. 4, we can see that the distribution is shifted to 
the right in case the distribution is non-thermalized. 
A partially non-thermalized strange quark distribu¬ 
tion will have more strange quarks with sufficient 
energy to form cc pairs or J/ip. This would lead to 
a non-negligible value of A q ^ cE and AThis 
would infuse more cc and J/ij} into the system, lead¬ 
ing to a slightly different value of suppression. Com¬ 
paring the experimental charm quark distribution 
and the modified Fermi Dirac distribution, there is 
a substantial difference in the population of charm 
quark at lower energy. Since the momentum range 
at ALICE forward rapidity and CMS mid rapidity 
is on the higher side, this difference should not im¬ 
pact the recombination calculation considerably for 
these two cases. The impact on ALICE mid ra¬ 
pidity could however be higher. It is also observed 
that shadowing calculation is a relatively uncertain 
quantity. Though we have not depicted the uncer¬ 
tainty in this work, EPS09 parametrization gives a 
range of uncertainty in the shadowing. 


Another source of uncertainty may be the for¬ 
mation of J/i[} during hadronization phase. To an¬ 
alyze this uncertainty, let us look at the formation 
of open charm mesons during the QGP phase and 
hadronization phase. The dissociation time of open 
charm mesons at the QGP deconfinement temper¬ 
atures is expected to be very small (less than a 
fraction of a fermi) compared to that of hidden 
charm mesons, and thus any open charm meson 
in the QGP medium will disappear quickly [28j. 
One may safely assume that open charm mesons 
would be formed mainly during hadronization. Dur¬ 
ing hadronization, due to the huge abundance of 
light quarks (possibly 2 or 3 orders of magnitude 
more than charm quarks), formation of open charm 
meson would be predominantly much more than 
the formation of J/-0 or any other hidden charm 
mesons. This can be seen from comparison of J/i/j 
and D meson production from pp collision exper¬ 


iments too @, 29], after adjusting for luminosity. 
As a result, the number of charm quark or anti¬ 
quark (=N C ) available for recombination to form 
J/ij) would reduce substantially. If one were to con¬ 
tinue with the quadratic dependence on available N c 
for recombination during hadronization, then the 
number of J/ijj formed due to recombination during 
hadronization would be negligible compared to the 
number of J/-0 formed due to recombination during 
the QGP phase. In this proposed model, we thus 
ignore the formation of J/'ijj during hadronization, 
which may lead to some discrepancy. 

Despite the above uncertainties, we see that our 
model has been able to explain the trend of central¬ 
ity dependent J/if) suppression data in the forward 
rapidity region obtained from the ALICE experi¬ 
ment and mid rapidity region from the CMS exper¬ 
iments with reasonable agreement. 

Many of the above mentioned discrepancies and 
uncertainties result in an error in either the initial 
conditions or in the dissociation or recombination 
rates. In order to estimate the consequence of esti¬ 
mation errors in the initial conditions and the dis¬ 
sociation or recombination rates, the J/^j suppres¬ 
sion is calculated by modifying the initial conditions 
and T diss and T recomb by 10%. In order to maximize 
the effect of error in initial conditions, the values of 
A r j/q(0) and IV c (0) are modified in opposite direc¬ 
tions, i.e. one is decreased by 10% and the other 
is increased by 10% and vice-versa. Similarly, on 
the dissociation and recombination rates, one is de- 














15 


creased by 10%, while the other is increased by 10%. 
Figure [L4] depicts the changes in final J/tp suppres¬ 
sion due to various modifications in the initial con¬ 
ditions, T diss and T recomb . It is interesting to note 



FIG. 14: Error propagation: The red curves 
indicate 10% change in initial conditions, while the 
blue curves indicate 10% change in 
dissociation/recombination constants. Detailed 
description of the curves is as follows: 
solid black curve: The calculated J/ip suppression; 
solid red curve: 10% decrease in Nj /^( 0 ) and 10% 
increase in 1V C (0); 

dashed red curve: 10% increase in Nj /^( 0 ) and 10% 
decrease in 1V C (0); 

solid blue curve: 10% decrease in Tdi SS and 10% 
increase in T recomb ; 

dashed blue curve: 10% increase in Tdi SS and 10% 
decrease in T recomb , 

that the final J/ip suppression is to some extent 
resilient to error in initial conditions and dissoci¬ 
ation or recombination rates. This can be under¬ 
stood by observing that the coupled rate equations 
lead to a negative feedback system. For example, 
if Nj/^( 0) were to increase by 10%, then a slightly 
larger number oi J/ip would dissociate into cc dur¬ 
ing the dynamical evolution and thus reducing the 
J/ip population. 

V. CONCLUSIONS 

We have calculated the recombination cross sec¬ 
tion at various temperatures for recombination of 


cc into J / ip at both ALICE and CMS experimen¬ 
tal conditions using pQCD. A temperature depen¬ 
dent Debye color screened phenomenological poten¬ 
tial has been utilized for computing the J/ip wave- 
function which is used in the calculation of the re¬ 
combination cross-section. We have established a 
set of rate equations which combines gluonic disso¬ 
ciation, collisional damping along with the recom¬ 
bination. CNM effects, namely shadowing is also 
incorporated in the current work. We find that 
the our final results explain the experimental J/ip 
suppression data at ALICE at forward rapidity and 
CMS at mid rapidity region to a good extent. This 
indicates that pQCD can be a reasonable approx¬ 
imation for recombination calculations at ALICE 
forward rapidity and CMS mid rapidity. At ALICE 
mid rapidity, where the momentum ranges are much 
lower, the discrepancy due to pQCD calculations is 
larger. However, pQCD is able to capture the trend 
that recombination would be higher at ALICE mid 
rapidity than at forward rapidity. A more exact 
estimate of recombination at ALICE mid rapidity 
would require non-perturbative QCD. We have also 
shown that for charmonium, modification of the for¬ 
mation time due to temperature leads to negligi¬ 
ble suppression due to color screening. The other 
causes of the difference between the experimental 
data and the predicted suppression could possibly 
be attributed due to the limitations in the accu¬ 
racy of determining the initial conditions for the 
rate equations, incomplete thermalization of strange 
quarks, limitation in the precise estimate of shad¬ 
owing effect and J/ip yield during hadronization. 
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Appendix A: Feynman Diagrams, Group 
Theoretic Factors and r recom b 

The Feynman diagrams pertaining to the calcu¬ 
lation of cjrecomb is shown in Fig. [[5l 







FIG. 15: cc scattering to emit a gluon and J/ij). 
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The group theoretic color factors appearing in 
the calculation of the above diagrams is given be¬ 
low. In the equations written below, G hJ refers to 
the group theoretic color factor appearing in the 
calculation of M* Mj. N c = number of colors = 3. 
the value of h is given by h abc = d abc + if abc where f 
is the totally antisymmetric structure constants of 
SU(3), while d is totally symmetric. 

G 1(1 = N C /A(N 2 - 1); 6 1)2 = iN c / 8 (N 2 - 1); 

Gi ,3 = iN c / 8(1 - N 2 C )-G 1A = iN c / 8 (N 2 - 1); 

G 1)5 = iN c /8(N 2 c + 1 ); Gb , 6 = i / 8 (N 2 + 1 ); 

Gg 7 = i/ 8 (N 2 - 1); G 1j8 = i / 8 (N 2 + 1); 

Gi 9 = i/ 8 (N 2 - 1); G l io = 2i.f ba ' c 'f bac f a ' ae h c ' ce + 
N 2 / 8 (l-N 2 ); 

G 2 , 2 = N*/ 8 -N c /A+l/ 8 N c -G 2>3 = 1/8(7V C +1/7V C ); 
G 2 4 = 4 h a ' ba h a ' ab \ G 2 5 = 4 h ba ' a h ba ' a ] 

6 2 ’ 6 = N 2 /8 + 1/4 1/8 N 2 C - G 2 ,7 = 1/8(1 - 1/N 2 ); 

G 2 , 8 = 1/8(1 - 1/N 2 )] G 2>9 = 1/8 (N 2 - 1/N 2 )] 

G 2 , 10 = i/ 8 (N 2 — 1 ); 

6 3 . 3 = 1V//8 - lV c /4 + 1 / 8 N C ] G 3 ,4 = Ah aa ' b h aa ' b \ 
G 35 = Ah ba ' a h baa 'n ; 

6 3 ’ 6 = 1/8(1 - 1/1V C 2 ); G 3 , 7 = 1V//8 + 1/4 - 1/8N 2 ] 
63 ,s = 1/8(1V C 2 - 1/1V C 2 ); 6 * 3,9 = 1/8(1 - 1/-/V 2 ); 

G 3 , 10 = i/ 8 (N 2 — 1); 

6 4 . 4 = Nl/ 8 -N c /A+l/ 8 N c -G^ = 1/8(N C + 1/N C ); 
G 4 , 6 = 1/8(1 - 1/1V C 2 ); 6 * 4,7 = 1V c 2 /8 + 1/4 - 1/81V C 2 ; 

6 4 . 8 = N 2 /8 + 1/4 - l/ 8 iV c 2 ; G 4 , 9 = 1/8(1 - l/N 2 ); 

64,10 — i/8(A^ 2 + 1); 

65.5 = iV c 3 /8-iV c /4+l/8iV c ; G 5)6 = l/ 8 (N 2 -l/N 2 )-, 
65,7 = 1/8(1 - 1/-/V 2 ); 6 5 , 8 = 1/8(1 - 1/N 2 )] 

65.9 = N 2 / 8 + 1/4 - 1/8 N 2 ; 65, 10 = i /8 (N 2 + 1); 
6 66 = TV 3 / 8 - N c /A + l/ 8 iV c ; 6 6 7 = Ah a ' ab h a ' ba \ 
6 6 ’ 8 = l/ 8 (iV c + 1 /N c ); 6 6 , 9 = Ah ba ' a h ba ' a - 

66. 10 = iN c / 8 (N 2 — 1); 

677 = TV 3 / 8 - N c /A + 1/8 N c ; 67 8 = Ah aa ' b h aa ' b ] 

67.9 = 1 / 8 ( 7 V C + l/AT c ); 67,10 = iN c /8(N 2 c + 1 ); 

Gs 8 = 7V 3 /8 - 1V c /4 + l/ 8 iV c ; 6 8 9 = 4/i a,f,a // ,a ' a ; 

6 8 .1 0 = * N C / 8 (N 2 — 1); 

69,9 = N*/ 8 -N c /A+l/ 8 N c - 6 9 , 10 = i N C / 8 (N 2 1); 

610.10 = N c / A ( N 2 - 1 ); 
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Appendix B: Infrared divergence of Fermion 
propagator 

This section proves the cancellation of the 
infrared divergence due to the fermion propagator. 
When the emitted gluon momentum approaches 0, 
the J /"0 center of mass reference frame becomes 
equivalent to the cc center of mass reference frame. 
The two momenta p\ and P 2 becomes equal and 
opposite. With this condition, let’s consider the 
sum of the Feynman diagrams 6 and 7. 


The term - r J — 

—m 


_ i( g 

k'f—m 2 2 p x .kg 


2 Px-kg 


. The first 


term will have no infrared divergence as k g —> 0 , and 
only the second term needs to be taken care of. Here 
p x is either p\ for the expression for Ad$, or p 2 in the 
case of M 7 . 


With the first term removed, the non scalar 
terms in the expression consist of 


{v(p2)(igiut a )v(k , ) S j 


= 


v{p 2 ){igiX)v(k')(- 


-w 


k\ 


u(k)(ig r ]f a t a 


kj i — m 




dke^ b (k g )( Bl) 


(u(k)(igYt a ) ^ 2 ^ k ( i 9d^ b ) u (Pi)^ j + 
( (v(P2){igint b )(ig'yX)v(k')^ 


(B 6 ) 


(B7) 


m 7= I (v(p2)(ign/vt b ) ^ m (ig^ u t a )v(k')^j 

(—p—) ( (u(k)(ig-fX)u(piX T (k, k')) ) 

dk e^ b (k g ){ B2) 


After contracting the o index, the Ad$ + M 7 
integrand would contain the expression 

((v(p 2 )(ig~/X)v(k ')) 

(u{k){ig'ft a ) -j^ — {igi»t b )u{pi)^j j (B3) 

+ (B4) 

( [v{P2)iigi, J t b ) ^ l _ m (ig^t a )v(k')^ 


From straightforward Dirac algebra, this is simpli¬ 
fied to 


f (y(p 2 ){ig'y v t a )v{tf) > ) 

(u{k)(ig'ft a )^^-{ig^)u{p 1 ) > j j + (B 8 ) 

fv(p 2 ) ( igt b ) (ig'yX)v(k') 

{u(k)(igYt a )u(p l ))'\ (B9) 


To see the infrared behavior, we write the 
particle spinor in terms of the antiparticle spinor 
and vice versa in expression for M 6 i.e. the first 
expression. 


(<‘(*0(ot 1 'O“(pi)) 


(B5) 


The first term is from the expression for M 6 , 
while the second is from the M 7 expression. The 
gluon polarization term is reintroduced in the end. 


v(k') = — 
v(p 2 ) = ~iu*(p 2 ) 72 
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We get 

v(P2) r Tuv(k') (BIO) 

= -iu*{p 2 )l 2 lv{-i)l 2 U*(k') 

= -(-l y* 2 u\p 2 ) lv u*{k') 

= -(-i)^V(p 2 )(-i)" =2 7 X(fc') 

= (u(p 2 )^u(k'))* 

(BH) 


The various minus signs used in arriving at Bll are 
as follows 

(a) — l 1 ^ 2 : from the commutation relation of 7 ^ and 
72 

(b) — : from 7 ! = — 1 

(c) — l u=2 : from 72 = —72 


space. Rotation of 7 r radians along a suitable axis 
interchanges p\ and p 2 as well as k and k'. 


= [(fi(fcb,«(p.)) (sMfrX*'))] (Bis) 

This expression for Mq is now identical to the 
expression for M 7 except for the scalar terms 2 ^ k 
and ^ 777 - Adding M e and M 7 and reintroducing 
the gluon polarization term, we get, 

()( p 1 ^(k,) + P 2 ^(k s ) 

Pl-kg p 2 • kg 

Since we are considering only transverse polariza¬ 
tion of the gluon, the time component is 0 and not¬ 
ing that pi = -p 2 , we get, 


On similar lines 

u(k)(Y 

2 pig , 


2 pi M 


u(pi 


2 p 1 .k g 


2 p\.kg 

v(k){Y)v(pi)) 


(B12) 


Taking the product of the two expressions Bll and 
B12, we get: 


= (afob„«(fc'))‘ (»(*0(7>(Pi))' (B13) 

This gives the expression for incoming particles 
of momenta k and k! and outgoing particles of ve¬ 
locity pi and p 2 . The conjugate transpose relates 
the above expression to the amplitude M for incom¬ 
ing particles of momenta p\ and p 2 and outgoing 
particles of momenta k and k' . Hence the above 
expression becomes, 


7 MPu^ikg) Pue i (k g ) 

(-)(——r- —r—) (B17) 

Pl-Kg P 2 -tig 

with i = 1,2,3. 

In the center of mass frame, if the gluon 3 momen¬ 
tum k g is perpendicular to the 3 momentum pj and 
p 2 , then pi.kg = p 2 .k g = E c E g , where E c = charm 
quark energy and E g = gluon energy in the center of 
mass frame. With this, the above expression evalu¬ 
ates to 0. When the gluon 3 momentum is parallel 
to p\ , e[kg) would be perpendicular to the 3 mo¬ 
mentum pi, and the numerator vanishes. Thus, if 
the gluon 3 momentum is perpendicular or paral¬ 
lel to charm quark (anti-quark) 3 momentum, the 
infrared divergence cancels. In other inclinations, 
there would be partial cancellation of varying de¬ 
grees. 

A similar analysis can be done for diagrams 2 
and 3, 4 and 5 and finally 8 and 9. 


= {u{k')^ u u{p 2 )) t {v(pi)(Y)v(k)y 

Zpi.rCg 

= (m(* 07„“(P2)) (®(Pi)(7*>W) 

= 4^7 Mk'hMpi)) (®(pi)(7>(*0)] (B14) 


The term in the square brackets is an inner prod¬ 
uct which would be an invariant to rotation in 3D 



